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Motivated by the statistical fluctuation of Dirac spectrum of QCD-like theories subjected to 
(pseudo)reality-violating perturbations and in the e-regime, we compute the smallest eigenvalue 
distribution and the level spacing distribution of chiral and non-chiral parametric random matrix 
ensembles of Dyson-Mehta-Pandey type. To this end we employ the Nystrom-type method to 
numerically evaluate the Fredholm Pfafhan of the integral kernel for the chG(0,S)E-chGUE and 
G(0,S)E-GUE crossover. We confirm the validity and universality of our results by comparing 
them with several lattice models, namely fundamental and adjoint staggered Dirac spectra of SU(2) 
quenched lattice gauge theory under the twisted boundary condition (imaginary chemical poten- 
tial) or perturbed by phase noise. Both in the zero-virtuality region and in the spectral bulk, 
excellent one-parameter fitting is achieved already on a small 4 4 lattice. Anticipated scaling of the 
fitting parameter with the twisting phase, mean level spacing, and the system size allows for precise 
determination of the pion decay (diffusion) constant F in the low-energy effective Lagrangian. 
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I. INTRODUCTION 

The understanding of the phase structure of fermion/gauge systems has posed a challenge for particle and nuclear 
physics communities alike. One possible resort to avoid the sign problem that the Euclidean Dirac operator becomes 
non-Hcrmitian and the Boltzmann weight complex in the presence of chemical potential /i is to substitute QCD by 
its hypothetical two-color version on a lattice [T], in which the Boltzmann weight (with pairs of degenerated quark 
flavors) is positive-definite even at finite fi due to the (pseudo)reality of the representations of SU(2). It has long 
been appreciated that the theories with quarks in the (pseudo)real representation exhibit exotic types of spontaneous 
breakdown of global flavor symmetry [2J|3]. Quarks and charge-conjugated antiquarks are combined into an extended 
Nambu multiplet, which in turn is expected to break down to the extended vector subgroup [2J. Then the effect of 
the chemical potential that breaks the extended flavor symmetry is unambiguously incorporated in the low-energy 
effective description through the flavor-covariant derivative [SHE]. The assertions and/or analytic predictions, possibly 
based upon the effective theory, can be quantitatively compared with the Monte Carlo simulations of SU(2) lattice 
gauge theory [9U12). provided that the lattice regularization respects the relevant flavor symmetry group. Accordingly 
the two-color QCD has served as a insightful testing ground for the realistic chromodynamics, as well as a tractable 
lattice model that is interesting by its own right. 

The difference in the global symmetries is a reflection of the difference in the antiunitary symmetries of the Dirac 
operators: Dirac operators in the fundamental representation of SU(2) and in the adjoint of SU(A^) are essentially real 
symmetric and quaternion selfdual, respectively, whereas that in the fundamental of SU(A?~ > 3) is merely complex 
hermitian, as Verbaarschot [T3] dubbed "threefold way" after Dyson's original proposal 14 . Thus the inclusion of 
the reality- or sclfduality- violating chemical potential /i in two-color QCD casts itself in statistical properties of Dirac 
spectra. As the symmetry-violating effect of /i in the two-color QCD is inherent in its low-energy effective theory and 
is well under control, one can predict the fluctuation of the Dirac eigenvalues in the e-regime (i.e. below the Thouless 
energy) from its zero-momentum part. This in turn is equivalent to the chiral Gaussian Orthogonal or Symplectic 
Ensemble (chGOE, chGSE) [T^J HB] in its non-Hermitian extended form by the introduction of schematic (real) \i 
component [HI [18]. 

At this point we should note that, the reality or sclfduality of the SU(2)-chromodynamic Dirac operator could as 
well be violated by the inclusion of any Hermitian component in the complex representation, most simply by U(l) 
electrodynamics or random phases, or even by a fixed Abelian Aharonov-Bohm flux background (i.e. twisted boundary 
condition or imaginary chemical potential) [191 - El] . These cases are distinct from the previously mentioned case in 
that the Dirac eigenvalues stay real even after the inclusion of symmetry violations and their statistical behavior 
exhibits crossover, rather than develop into the complex plane. In terms of the effective cr-model description these 
two cases are almost identical, save for the difference of the sign of tr BQ^BQ term (/x 2 F 2 or (ifi) 2 F 2 , see sect. IV). 

Crossover between universality classes of Hermitian random matrix ensembles [22TI24) . namely GOE-GUE and 
GSE-GUE, is extensively studied in the context of disordered [551 and quantum-chaotic Hamiltonians [2"Trf2l)] 
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with its time-reversal invariance slightly broken by weak magnetic field or AB flux applied [30]. On the other hand, 
the chiral or superconducting variant of universality crossover appears to be a relatively unexplored field so far. 
Previous attempts in this area either focused on the level number variance and the spectral form factor (both of 
which are integral transforms of the two-level correlator and insensitive to chiralness) of two-color QCD with AB 
fluxes versus GOE-GUE crossover [T!5], or have fruited in a series of tours-de- force by Damgaard and collaborators 
[3T] [32] devoted to the analytical computation of the level density and individual small eigenvalue distributions for the 
spectral crossover within the chiral Gaussian Unitary Ensemble (chGUE) class, due to the imaginary isospin chemical 
potential. Despite that analytic results for the microscopic spectral correlation functions are known for some time 
for chGOE-chGUE and chGSE-chGUE crossover [33] [31], they are yet to encounter with physical application. To 
the best of our knowledge, the only example of the crossover involving different Hermitian chiral universality classes 
discussed in a physical setting is the CI-C transition for the normal/supcr-conducting hybrid interface in a magnetic 
field |35j . Part of the aim of this paper is to present novel examples of physical application of the crossover between 
chiral Hermitian universality classes in the realm of lattice gauge theory, despite within rather academic settings of 
QCD-like theories. 

In this paper we shall show that such Dirac spectra indeed exhibit symmetry crossover from the chGOE or chGSE 
to the chGUE universality class, precisely as predicted by the chiral variants of parametric random matrix ensembles 
of Dyson [22] and Mehta-Pandey [23j [24] , both in the zero-virtuality region and in the spectral bulk. Through the 
excellent one-parameter fit to the parametric random matrix results, we extract the pion decay (diffusion) constant 
in the effective Lagrangian. As our method adopts the level spacing and the smallest eigenvalue distributions that 
are extremely sensitive to the fitting parameter as primary fitting observables (see Refs.[36 38 for recent efforts along 
this line), it enjoys a clear advantage over the methods using n-level correlation functions, and presents promising 
applications in analyzing the numerical data of QCD-like theories. 

This paper is composed as follows: In section II we briefly review the universality crossover in the spectrum 
of random matrices. Then we shall compute and plot the level spacing distribution in the bulk and the smallest 
eigenvalue distribution at the hard edge (origin) of the spectrum using the Nystrom-type method, the latter being our 
new contribution. In section III we measure the level spacing and smallest eigenvalue distributions of fundamental and 
adjoint staggered Dirac operators of SU(2) quenched lattice gauge theory with weak AB flux or phase noise included, 
and fit the spectral data with the predictions of (ch)GSE-(ch)GUE and (ch)GOE-(ch)GUE crossover. We show that 
an excellent one-parameter fitting can be achieved for all of our cases of concern, so that an accurate determination of 
the crossover parameter is possible. In section IV we extract the pion decay constant in the effective chiral Lagrangian 
from the flux dependence of the crossover parameter. We shall conclude in section V with some discussions on the 
two-color QCD+QED simulation and on possible directions of future study. 

II. PARAMETRIC (CHIRAL) RANDOM MATRICES 

Transition of spectral fluctuation from one universality class to another with a different antiunitary symmetry 
has been proved to occur both in quantum-chaotic and disordered Hamiltonians with weakly broken time-reversal 
invariance. The result is universal in a sense that the local spectral fluctuation is sensitive only to a single crossover 
parameter p defined below. The reason for this universality is traced back to the nonlinear a model governing the 
spectral statistics, which can either be derived by the conventional disorder averaging of random Hamiltonians |25j 
or by the summation over Sieber-Richter-like encountering multiplet of periodic orbits of chaotic dynamical systems 
[25], completely irrespective of the details of dynamics. Accordingly one can resort to the simplest model which yields 
the identical a model, i.e. parametric random matrix ensembles. Below we collect established results on the spectral 
correlation of the parametric random matrices of non-chiral and chiral types for completeness, and refer the reader to 
the original references [22 -24, 33 for their derivations. 

A. GOE-GUE and GSE-GUE crossover 

We consider an ensemble of N x TV Hermitian complex (quaternion) matrices H = Hg + iaH^, with H§ real 
symmetric (quaternion selfdual) and Ha real antisymmetric (quaternion antiselfdual) matrices distributed according 
to Gaussian measures of variance a 2 . This parametric (also called Brownian-motion or dynamical) random matrix 
ensemble interpolates between the two limiting cases, GOE (GSE) at a = and GUE at a = 1. Take a point A from 
the bulk part of the spectrum of H and denote the mean level spacing around A by A(X). Then the n-point correlation 
function of the eigenvalues {A;} of H in the vicinity of A is given, in the limit N — > oo, a — > and p = aa j A(X) fixed 
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(we follow Mehta's book 24 for the definition of p), as a Pfaffian, 

R n (%i, ■ ■ ■ ,x n ) = Pf(Z [K (xi,Xj)]" j=1 ) = .y/det [K(xi,Xj)]" j=1 , 
S(x,y) I(x,y) 



K(x,y) 



, Z = ia 2 ® In, (1) 



D(x,y) S(y,x) 

where Xi = Xi/A(X) are the unfolded eigenvalues, and S(x, y), D(x, y), I(x, y) as functions of r = x — y are given by 

[GOE-GUE] S(r) = , D(r) = - f duti^ sinur , I(r) = - f°° ^V 2 ^ 2 sinw , (2) 

nr it J irj, t) 

CiTj 7TT* 1 2 2 1 f\ r \) 2 2 

[GSE-GUE] S(r) = , D(r) = dvve~ 2p v sinvr , I(r) = - — e 2p v sinvr . (3) 

nr 7T J n 7T Jo v 

The probability -E(s) that an interval of width s contains no eigenvalue is then given as the Fredholm Pfaffian, or 
square root of the Fredholm determinant [25] : 



(-1)" 



E(s) = ^ ^—^ I dXf- I dx n Rnixx, ...x n ) = ^/Det(l - K s ) (4) 

n=0 



where K s is an integral operator of convoluting with the "dynamical" sine kernel K(x,y) (jf), fah, fcfy, restricted to 
the interval [0, s]. The probability distribution P(s) of level spacings s = Xi + \ — xi is given by its second derivative 
P(s)=E"(s). 

We should emphasize that for parametric random matrix ensembles, local correlations of unfolded eigenvalues in 
the vicinity of A still depend upon the mean level spacing A(X) of the eigenvalue window in concern through the 
parameter p. Accordingly, if the parameter a is adiabatically increased from zero to unity, universality crossover from 
GOE or GSE to GUE takes place at a different rate in each window in the spectrum (the denser the eigenvalues, 
the faster the speed of crossover). Also note that the universal intermediate behavior of spectral fluctuations appears 
only in a double limit where the system size tends to infinity and the a parameter to zero in a correlated manner; a 
simple thermodynamic limit N oo with a fixed would drive the whole spectrum to the GUE class. 

B. chGOE-chGUE and chGSE-chGUE crossover 

The chiral version of the parametric random matrix ensembles are simply obtained by setting the N/2 x A^/2 
block-diagonal parts of H = Hg + iaH^ to zero. Accordingly the matrix in concern takes the form 

rr ( H X + iaH 2 \rrrjNN 

H = I ,jj ictiy ) T - D J ' 1 '-" 2 : ~2~ x ~2~ ( c l ua t ernlon -) rea l matrices, (5) 

distributed according to Gaussian measures of variance a 2 . This ensemble interpolates between the two limiting cases, 
chGOE (chGSE) at a — and chGUE at a = 1. Since the nonzero eigenvalues of H occur in the ± pairs of equal 
magnitude, if suffices to retain only non-negative eigenvalues. The rt-point correlation function of the eigenvalues 
{Xi} of H in the vicinity of the origin is similarly expressed, in the limit N — > oo, a — > and p = aa/ A(0) fixed, 
as a Pfaffian Q with S, D, I given by (after substituting X —> tt 2 x 2 ,Y —> ir 2 y 2 into the "Laguerre-type" formulae 
(7.2.28~56) of Ref.[33] and multiplying them by 2ny/xy), 



[chGOE-chGUE] S(x,y) = Z M™)vM*v) + joM f°° dve -W-*>) Jq{vx) 

{ x 2 -y 2 2 J w 



D(x,y) = -^l I dvv 2 e 2pV {x.h(vx)J a {vy) ~ J (vx)yJ x {vy)} , 



DC 



I(x,y) = ^ J dvvj^ due-P v > {J (vux)J Q (vy) - J Q (vx)J (vuy)} , (6) 

[chGSE-chGUE] S {x,y)=nVxy( ^^ 

I x 2 -y 2 2 J 





x z 






dv v 


[ due" 




'0 


dv v 2 





D(x,y) = ^ 

L Jo 

I(x,y) = ^ [ dvv 2 e- 2 »~ v ~ {xJ x {vx)J {vy) - J {vx)yJi{vy)} , (7) 
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where Xi = Xi/A(Q) are the unfolded eigenvalues (see also [5i]V 

The probability E(s) that no eigenvalue is smaller than s is again given as the Fredholm determinant Q, where 
K s in this case is an integral operator of convoluting with the dynamical Bessel kernel K (x, y) |l]), (F|, restricted 
to the interval [0, s\. The probability distribution p±(s) of the unfolded smallest eigenvalue s — Ai/A[Q), which is one 
half of the very central level spacing, is then given by its first derivative pi(s) = —E'(s). 



C. Nystrom-type method 



A simple but exceptionally efficient way of evaluating the Fredholm determinant of a trace-class integral operator 
K s acting on the Hilbert space of L 2 -functions / over an interval [0,s] by (K s f)(x) — J* dy K (x , y) / (y) , such as the 
one with the sine kernel S(r), is the Nystrom-type method [351 SO]- It simply discretizes the Fredholm determinant: 

Det(l - K s ) ~ det [<% - K(xi,Xj)y/w^w]]™. =1 . (8) 

Here the quadrature rule {x\, . . . , x m ; w\, . . . , w m } consists of a set of points {xi} taken from the interval [0, s] and of 

PS m 

positive weights {wi} such that / f(x)dx ~ f(xi)wi. Efficient choices for the quadrature rule are Gauss (sampling 

J ° ' i=l 

at the Legendre nodes) and Clenshaw-Curtis (sampling at the Chebyshev nodes), for which the computational cost 
grows optimally as 0(m 2 ) and O(mlogTO), respectively [H]. As the order to of the approximation increases, the RHS 
of© is proven to uniformly converge to its LHS. The convergence is rapid and exponentially- fast; the approximation 
error decays as O(e~ const m ) [3D]. In this paper we choose the Gauss quadrature rule, as 15-digit accuracy is already 
attainable only with to = 5, for the Fredholm determinant P(0.1) for the sine kernel. An extension to the matrix- 
valued kernel ([I]) is trivial: one merely takes the determinant over the matrix indices as well. Practical significance 
of the method in the context of random matrices and stochastic processes is recently re-appreciated and stressed in 
Ref-gDJ. 

We have applied this Nystrom-type method to the dynamical sine kernel 0, (j2J), ^ and the dynamical Bessel 
kernel Q, ©, to obtain P(s) and pi(s) for G(0,S)E-GUE and chG(0,S]E-chGUE crossover, respectively. In 
order to achieve accuracy that is needed for computing the first or second derivatives (pi(s) or P(s)) to a good 
precision, we have chosen the approximation order m to be (at least) 20 for the former and 100 for the latter, and 
confirmed the stability of the results under the increment of to. Numerical results for the region < s < 3 ~ 4 and for 
the parameter range p < 1 are exhibited in Figs. 1 and 2 (left) (P(s) for non-chiral random matrices) and in Figs. 3 
and 4 (left) (pi(s) for chiral random matrices) 1 Although Mehta-Pandey (in their second paper of Ref. 23J) have 
expressed the Fredholm determinant for the dynamical sine kernel in terms of eigenvalues of an infinite-dimensional 
matrix (each matrix element of which is an integral involving prolate spheroidal functions), these numerical plots of 
P(s) and pi(s) for (namely the chiral version of) parametric random matrix ensembles do not seem to have appeared 
explicitly in the literature, to the best of our knowledge 2 . Also plotted in the figures are P(s)'s for the two limiting 
cases p = and p = oo, i.e. for GOE (GSE) and GUE obtained by Jimbo-Miwa-Mori-Sato [12] in terms of a solution 
to the Painleve V transcendental equation subjected to an appropriate boundary condition, and pi(s)'s for chGOE 
(chGSE) and chGUE [43], 

exp (— f s) exp i 

Pl (s) = -~J coshfs exp(-^fs 2 ) [chGSE] . (!)) 



d 



exp i 



— s 2 

8 * , 


| [chGOE] 


8 * , 


) [chGSE] 


4 J 


| [chGUE] 



We immediately observe the asymptotic behaviors of P(s) and pi(s) for parametric (chiral) random matrices at 
finite p, 

P^( S )~C p .s 2 , p^^-cps 1 ( S «l); logPW(s), logp ( f \s)~- lp s 2 (s > 1). (10) 



1 Due to Kramers degeneracy, a half of all the level spacings of GSE are zero and its distribution has a peak at the origin, Pgse(* ~ 0) ~ 
5(s)/2. Thus its smooth part is rescaled and normalized as f^° Pgse{s)<1s = 1/2. Accordingly P(s) for GSE-GUE crossover becomes 
peaky near s w as the p parameter is decreased. It also leads to the loss of accuracy of the Nystrom approximation at very small p. 
The first reference of 30 did not use the analytic form of P^(s) for GOE-GUE crossover plotted in Fig.l (left), but employed the 
Wigncr-surmised form (see also 1531 ). 
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FIG. 1: Level spacing distribution P(s) and two-level correlation function R,2(x,0) for GOE-GUE crossover. The p parameter 
ranges 0.02 < p < 0.70 or 0.50 by step 0.02. The two bounding curves correspond to the GOE and GUE limits, p — and oo. 





FIG. 2: Level spacing distribution P(s) and two-level correlation function R2(x,0) for GSE-GUE crossover. The p parameter 
ranges and 0.01 < p < 0.70 or 0.50 by step 0.01. The two bounding curves correspond to the GSE and GUE limits, p — and 
oo. 





FIG. 3: Smallest eigenvalue distribution pi(s) and microscopic level density Ri(x) for chGOE-chGUE crossover. The p 
parameter ranges 0.02 < p < 1.00 by step 0.02. The two bounding curves correspond to the chGOE and chGUE limits, p — 
and oo. 
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FIG. 4: Smallest eigenvalue distribution p\ (s) and microscopic level density R\ (x) for chGSE-chGUE crossover. The p parameter 
ranges 0.02 < p < 1.00 by step 0.02. The two bounding curves correspond to the chGSE and chGUE limits, p — and oo. 



Here C p ,c p ,7p are p-dependent constants, monotonically varying in p, in ranges 7r 2 /3 < C p < oo, 7r 2 /2 < c p < oo, 
tt 2 /16 < 7 P < tt 2 /8, i.e. in between the (ch)GUE and the (ch)G(0,S)E limits. 

For a bookkeeping purpose, we exhibit plots of the two-level correlation function i?2(a;, 0) for non-chiral cases in 
Figs. 1 and 2 (right), and the single-level density Ri(x) for chiral cases in Figs. 3 and 4 (right). Note that the first peaks 
of i?2(z,0) and Ri(x) (right figures) are comprised of the corresponding P(s) andpi(s) (left figures), respectively. In 
the subsequent section, these analytic results will be tested to fit the Dirac eigenvalue data numerically obtained from 
modified SU(2) lattice gauge models. The practical advantage of adopting distributions of individual level spacings 
(P(s) and Pi(s)) over n-level correlation functions (i?2(£,0) and Ri(x)) for fitting is clear from the figures. As the 
oscillation of the latter consists of overlapping of multiple peaks, the characteristic shape of each peak is inevitably 
smoothed out, leaving us with a rather structureless curve for which an accurate fit is difficult. On the other hand, 
the former and its cousins Pk(s) (the distribution of the k th smallest eigenvalue [4"5] ) 



Pk{s) 



d 
ds 



(fc-1)! 



d_ 

dz 



k-l 



Det(l - zK s ) 



(11) 



that could as well be computed by the Nystrom-type method are very sensitive to the value of the crossover parameter, 
because Eqs. (10) imply that the ratio of P{s) or pi(s) for the orthogonal and sympletic classes to that for the 
unitary class grows as exp E j^- for large s (it exceeds 10 at s ~ 2, 10 2 at s ~ 2.7, and 10 4 at s ~ 4). Therefore 

P^(s) and p±{s) should in principle admit very sharp one-parameter fitting by the least square method or simply 
from the tail of the curve (in the range of s where the systematic deviation due to finite size is not prominent) , as 
done for the Anderson tight-binding Hamiltonians at the metal-insulator transition 44 versus the critical random 
matrix ensembles interpolating G(0,U,S)E and Poisson statistics [I5J. In the chGSE-chGUE case, the microscopic 
level density Ri(x) (Fig. 4b) would also be suited for fitting, due to its characteristic oscillatory behavior for a wide 
range of x. 



III. DIRAC SPECTRUM 

A. Antiunitary symmetry 

Dirac operator for the "real" QCD, i.e. for quark fields belonging to the real or pseudoreal representation of 
the gauge group, is known to possess a particular antiunitary symmetry unlike that for the complex representation 
[15] . Namely, Euclidean Dirac operator D = 7 P («<9 A1 + A®T a ) for quarks in the fundamental representation of SU(2) 
commutes with Ct 2 K. Here C is the charge conjugation matrix satisfying Cj^C^ 1 — —7^, t 2 is one of the generators 
of the SU(2) gauge group satisfying T2T a T 2 ~ 1 = — r*, and K denotes the complex conjugation. As (CT2K) 2 = +1, D 
can be brought to a real symmetric matrix by a similarity transformation. Similarly, Dirac operator for quarks in the 
adjoint representation of SU(A?") commutes with CK. As (CK) 2 = —1, D can be brought to a quaternion selfdual 
matrix. By the same token, Dirac operators in fundamental representations of O(A^) and Sp(2iV) gauge groups are 
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real symmetric and quaternion selfdual, respectively. It is also well known that the reality and the selfduality of 
t he c ontinuum Dirac operators are interchanged for the corresponding Kogut-Susskind staggered Dirac operators 



((15) below), i.e. D KS in SU(2)-fundamental is essentially quaternion selfdual whereas D KS in SU(iV)- adjoint is real 
symmetric, due the the absence of the charge conjugation matrix |46j . The spectrum of SU(2)-fundamental staggered 
Dirac operator was indeed numerically shown to belong to the chGSE class [47] . 

Since the Dirac operator in the fundamental representation of U(l) (or SU(iV' > 3) if one prefers) in continuum or 
on a lattice possesses no such antiunitary symmetry, Dirac operator in SU(2)-fundamental (SU(-/V)-adjoint) xU(l), 

D = lll {id ll +AlT a + B li ), (12) 

where T a 's are the generators of the corresponding representation of SU(2) or SU(iV), and the U(l) gauge field, 
either weakly fluctuating or fixed as a background, is endowed with a weakly broken antiunitary symmetry as compared 
to the pure SU(2)-fundamental or SU(iV)-adjoint case. 

B. Modified SU(2) lattice gauge theory 

In contrast to the parametric random matrices H = H$ +iaHx for which the breaking of the antiunitary symmetry 
is uniquely parametrized by a — y (||Im_ff 1 | 2 )/ (||Rc H || 2 ) , it is not straightforward to identify the bare U(l) gauge 
coupling constant as proportional to a, due to the gauge invariance: a Dirac operator that appears complex could 
actually be a U(l) gauge transform of some purely real matrix. In that case its spectral fluctuation would perfectly 
be described by chGOE or chGSE. Thus in this paper we restrict ourselves to simpler models without such a subtlety: 
SU(2) quenched lattice gauge theory associated with twisted boundary condition or coupled with U(l) noise, each of 
which does contain an unambiguous counterpart of the a parameter. 

We consider SU(2) quenched lattice gauge theory under the twisted boundary condition (TBC), that is to multiply 
SU(2) links variable at the temporal boundary of the hypercubic lattice of size V — L A by a constant phase 3 

e* 2 ^ (n 4 = L, fi = 4) 

1 (else) ' (L6) 

with (p <C 1. This twisting is gauge-equivalent to the fixed vector U(l) background = (2mp / 'L)<5 M- 4 of AB flux 
[30] or the imaginary chemical potential /i = i2nip/L. The flux ip is the measure of antiunitary-symmetry violation and 
plays the role of a in the parametric random matrices. Its effect on the low-energy effective Lagrangian is completely 
dictated through the flavor-covariant derivative [20] . just as in the case of symmetry- violating real chemical potentials 
0. 

We also consider SU(2) + U(l) phase noise (PhN) model, where an i.i.d. random phase e l0n -^ taken from the uniform 
distribution 6 £ [— p,p] (p <C 1) (one could alternatively adopt Gaussian distribution of variance p 2 ) is multiplied 
to each link SU(2) variable. The parameter p of randomness is expected to be proportional to a in the parametric 
random matrices, since 

(||Imi?|| 2 ) _ (||[/sing|| 2 ) _ 1 

(||ReC|| 2 ) - (||t/cosc?|| 2 ) - K ' 3 P ' 1 ' 

The PhN model was previously used to uncover the fake nature [J5] of the apparent chGUE-chGSE crossover of 
the Dirac spectrum of Ginsparg- Wilson type |49| . so it would be a good exercise to exploit it again to our actual 
universality crossover. 

For our purpose of confirming that the lattice model exhibits chG(S,0)E-chGUE crossover, it is sufficient to focus 
on the strong coupling region of SU(2) at /3 = 4/ g 2 = ~ 1.5, where the smoothed level density at the origin, 1/A(0), 
is sufficiently above zero and the chiral symmetry is spontaneously broken, accepting that the model is away from the 
continuum limit. Accordingly we employ the simplest algorithm possible: the unimproved plaquette action and the 
10-hit heat-bath update coupled with overrelaxation. The staggered Dirac operator 

D KS — - (n S U e i0n ^-n S ^-U^ P ~ l9m ^) n = C-l )ni+-+n„-i /ir\ 



3 As the Dirac operator possesses (pseudo-)reality either for periodic or anti-periodic boundary condition, we adopt periodic conditions 
for all directions for simplicity and consider a small deviations from it. 



is diagonalized using the standard LA package. 

Because of our need to detect possibly small deviations of spectral fluctuation from the universal random matrix 
statistics at cither end (p — or oo), we give priority to the number of independent samples and perform our 
simulation on a lattice of the smallest size V — 4 4 . This choice is sufficient for measuring the short-distance behavior 
of eigenvalues within up to 3~4 mean level spacings and determining the p parameter precisely. In this region the 
systematic deviation due to the small size of the lattice is less prominent (they shall manifest at larger separation) 
than the statistical fluctuation. Only if one dares to test the symmetry crossover on the large-distance correlation of 
eigenvalues, such as the number variance S 2 (s), the use of larger lattices would become essential. 

C. Fitting Dirac spectra 

Our procedure of fitting Dirac spectral statistics to the parametric random matrix predictions consists of the 
following steps: 

1. Determination of -4(0). Perform the pure SU(2) simulation for each j3 and measure the fundamental and 
adjoint staggered Dirac spectrum {A^} for O(10 5 ) configurations. Taking for granted that the small Dirac 
eigenvalues obey the chG(S,0)E statistics, determine the mean level spacing at the origin A(0) by the following 
two ways: 

(a) Using the least square method, fit the histogram of the smallest Dirac eigenvalues pi(Ai) to the rescaled 
chG(S,0)E result p x {Xi/A)/A, by varying A. The fitting range is chosen to be within 34 (chGSE) and 
2.2A (chGOE), which is divided into bins of widths 0.14. Errors of A are estimated as those increase \ 2 by 1. 

(b) Compare the mean value of the smallest Dirac eigenvalue (Ai) with that of the (unfolded) chiral random 
matrix eigenvalues from 

{xi)= I s Pl (s)ds = J— erfc^ [chGOE] , J— [chGSE] =*> A(0) = (16) 

JO V 7T V2 V 7T (Xl) 

Typically the values of mean level spacing determined by the above two methods agree within 0.1%. This 
provides a ground for a precise determination of the crossover parameter at the origin in the step 2. 

2. Fitting the smallest eigenvalue. Multiply SU(2) link variables U n ^ by the phases e l6, "- fl (either TBC or 
PhN), and measure the Dirac spectrum {A;} for N con { = O(10 4 ) independent configurations. The unfolded 
smallest eigenvalue is still defined by x\ — Ai/Z\(0) with respect to A(Q) determined in the step 1, i.e. from 
unmodified SU(2) configurations. Fit the frequencies of Xi to the prediction by the least square method as 
follows: we choose the valid range of fitting the smallest eigenvalue X\ to be [0, 2.8] (fundamental) and [0, 
1.6] (adjoint), and divide them into r = 20 segments Iq, . . . ,I r -i of equal width. The reminder is denoted by 
I r = [2.8 or 1.6, oo). Define Xp = Hl=o — /» )V/i ) where Fi is the frequency of x\ falling on the i th 
interval It and — N ccm { Jj p[ p \s)ds its expected value from the chG(S,0)E-chGUE parametric random 
matrices. By varying p, find the optimal value of p that minimizes Xp an d estimate the error of p as that 
increases Xp by 1. 

3. Fitting level spacings. Fluctuation of Dirac eigenvalues in the spectral bulk is by itself not directly related to 
the low-energy effective description of the gauge theory, but rather to the diffusion of quarks in a hypothetical 
'time' evolution with Dirac operator as its Hamiltonian [50]. Nevertheless we shall argue that the distribution 
of level spacings from the spectral 'plateaux' [A m , Am] adjacent to but not including the origin, in which the 
mean level spacing A(X) is well approximated as a constant close to 4(0) (Fig. 5), provides an efficient way of 
determining F. It also serves as a clearer indicator of the presence of the universality crossover in modified SU(2) 
lattice gauge theories and as a consistency check of the p value directly measured from the smallest eigenvalues 
distribution. In order to avoid possible distortion of the level spacing distribution due to the chiralness, we take 
A m not too close to the origin and set it to be, at smallest, the 11 th eigenvalue. 

Define unfolded level spacings in that eigenvalue window as s = (Xi+i — Xi)/A(X), Ai,Ai+i € [A m ,Aiv[]- In case 
that the density profile near the origin is not quite flat (such as D in SU(2)-fundamental at (3 = 0.5 ~ 1.5), 
we alternatively used unfolding by the linearized level density p(X), A^ — > x, L = J p(X)dX even within the 
window chosen, to offset the explicit local variation in the mean level spacing. For a window containing ~ 20 
eigenvalues, the distortion in Pdata(s) caused by the small local inhomogeneity of the p parameter through its 
dependence on A = 1/p is expected to be negligible. Find the best fit of the histogram of the unfolded smallest 
Dirac eigenvalues -Pdata(s) to the crossover chiral random matrix result P^(s) (Figs. 2a or la) by varying p, as 
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FIG. 5: Macroscopic level densities p(A) of SU(2)-fundamental (left) and adjoint (right) D at various (3 in the strong coupling 
region. Lattice size: 4 4 , number of configurations: 100000. Level spacings from the eigenvalue windows (plateaux) [A m , Am] 
marked by orange strips in the vicinity of A = are used for fitting P^ p \s). 



done in the step 2. Thanks to the enormous gain in statistics due to the spectral averaging, we can safely set 
the fitting range to be as large as [0, 4] and divide it into r = 40 segments. 



D. Simulation results 



We present the outcome of one-parameter fitting of smallest eigenvalue distributions (SED) and level spacing 
distributions (LSD) to the parametric random matrix ensembles. We set the symmetry-breaking parameters to be 
if = 0.01 ~ .06 for the SU(2)+TBC model and p = 0.02 ~ .10 for the SU(2)+PhN model, and performed simulations 
at the gauge coupling constants (3 = 0, 0.5, 1, 1.5 on a lattice of dimensions V = 4 4 . N con f = 40000 configurations are 
generated for each set of parameters. Optimal values of the p parameter for the TBC (PhN) model in fundamental 
(F) and adjoint (A) representation are tabulated in the six (seven) columns in the middle of Table I (II) . 

Sample plots of pi(s), Ri(x), and P(s) for /3 = are exhibited in Figs.6-8 (SU(2)+TBC) and in Figs.9~ll 
(SU(2)+PhN). In all figures, measured histograms are plotted by colored dots, and parametric random matrix results 
that fit the data optimally are shown in curves of the same color (p = .01 or p = .02 in red to p = .06 or p = .10 
in green). Also plotted in the figures are the results from chG(S,0)E or G(S,0)E (black real line), and chGUE or 
GUE (broken line). Note that the microscopic level density Ri(x) is not fitted to the corresponding data by the least 
square method, but is merely substituted the p parameter determined from the SED. 

Even at inspection, one is convinced the accuracy of one-parameter fitting in all cases. The x 2 /d.o.f. deviations 
from optimally-fitting random matrix distributions are summarized in Table III. Considering the smallness of our 
lattice (4 4 ), the goodness-of-fittings achieved is astonishing (except for LSDs for SU(2)-fundamental at very small 
p < 0.1, where the distribution becomes extremely peaky at small s due to the onset of Kramers degeneracy and the 
fitting error is inevitably enhanced). These listed values are comparable to those reported in the pioneering papers 
in this field [31] using chGUE-chGUE (i.e. Hermitian) crossover, fitted to the spectral data from larger lattices: 
X 2 /d.o.f.= 0.33 for quenched QCD on 12 4 and x 2 /d.o.f.= 1.13 ~ 1.33 for dynamical QCD on 6 4 . 

From Tables I and II one immediately notices that for a fixed bare coupling (3, tp-p and p-p plots arc all quite 
linear. In Figs. 12 we exhibit sample plots of VAp/tp for SU(2)+TBC and y/Ap/p for SU(2)+PhN 4 , showing that 

(i) these ratios, namely from LSDs, are very stable under the change of symmetry-violating parameters ip or p, and 

(ii) the values of p determined from SED and from LSD are in a good agreement, as should be. This precisely 
linear dependence is essential in accurately determining the pion decay constant from SU(2)+TBC models in the next 
section. We note that in the strong coupling limit (3 = 0, four plots nearly overlap on the top of each other, regardless 
of the representation being fundamental (F) or adjoint (A) (also see Figs. 13). The splitting between fundamental and 
adjoint becomes more apparent at larger value of (3. 



4 The factor \f~A is included to facilitate extraction of F 2 /H in the next section. 
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TABLE I: SU(2) + TBC model on V = 4 4 : Crossover parameters and low-energy constants 



p 


rep/dist 




A 


E 


P 


/~a i 


F 2 /E 


T-l2 

F 






r \ \ l 

[A m , AmJ 






9?=.01 


.02 


.03 


.04 


.05 


.06 











F/SED 




.00933(1) 


1.315(2) 


.061(1) 


.123(1) 


.181(1) 


.245(2) 


.305(2) 


.364(2) 


.375(1) 


.221(1) 


.290(2) 




F/LSD 


[.10, .30] 


.00929(0) 


1.321(0) 


.061(0) 


.122(0) 


.183(1) 


.243(1) 


.303(2) 


.364(4) 


.374(1) 


.219(1) 


.290(1) 




A/SED 




.00617(0) 


1.989(1) 


.076(1) 


.153(2) 


.224(3) 


.309(5) 


.391(8) 


.453(11) 


.382(2) 


.229(3) 


.455(6) 




A/LSD 


[.10, .50] 


.00618(0) 


1.984(0) 


.075(2) 


.150(4) 


.224(6) 


.297(8) 


.371(14) 


.447(25) 


.374(5) 


.219(5) 


.435(11) 


0.5 


F/SED 




.01020(1) 


1.203(2) 


.055(1) 


.110(1) 


.160(1) 


.210(2) 


.263(2) 


.313(2) 


.339(1) 


.180(1) 


.217(1) 




F/LSD 


[.10, .26] 


.01004(0) 


1.223(0) 


.052(0) 


.103(0) 


.155(1) 


.206(1) 


.257(2) 


.306(2) 


.329(1) 


.170(1) 


.208(1) 




A/SED 




.00626(2) 


1.960(6) 


.067(1) 


.131(2) 


.200(3) 


.263(4) 


.327(6) 


.386(8) 


.332(2) 


.173(2) 


.339(5) 




A/LSD 


[.10, .50] 


.00622(0) 


1.972(0) 


.067(3) 


.134(4) 


.202(5) 


.269(7) 


.333(11) 


.407(17) 


.337(4) 


.179(5) 


.353(9) 


1.0 


F/SED 




.01150(5) 


1.067(5) 


.056(1) 


.102(1) 


.148(1) 


.193(1) 


.241(2) 


.286(2) 


.329(1) 


.170(1) 


.181(1) 




F/LSD 


[.10, .30] 


.01105(0) 


1.110(0) 


.048(0) 


.095(0) 


.143(1) 


.190(1) 


.237(1) 


.284(2) 


.319(1) 


.160(1) 


.177(1) 




A/SED 




.00631(3) 


1.944(9) 


.065(1) 


.132(2) 


.200(3) 


.266(4) 


.338(6) 


.397(8) 


.336(2) 


.177(2) 


.344(5) 




A/LSD 


[.10, .50] 


.00632(0) 


1.943(0) 


.066(3) 


.133(4) 


.199(5) 


.266(7) 


.331(10) 


.394(16) 


.335(4) 


.176(4) 


.343(9) 


1.5 


F/SED 




.01412(1) 


0.869(0) 


.061(1) 


.101(2) 


.138(1) 


.178(1) 


.218(2) 


.257(2) 


.329(1) 


.170(1) 


.148(1) 




F/LSD 


[.13, .28] 


.01280(0) 


0.959(0) 


.043(0) 


.085(0) 


.127(1) 


.169(1) 


.210(1) 


.251(1) 


.319(1) 


.160(1) 


.153(1) 




A/SED 




.00654(1) 


1.877(3) 


.067(1) 


.128(2) 


.193(3) 


.251(4) 


.312(6) 


.392(8) 


.336(2) 


.177(2) 


.333(5) 




A/LSD 


[.10, .50] 


.00650(0) 


1.887(0) 


.065(3) 


.128(4) 


.194(5) 


.257(7) 


.322(10) 


.386(15) 


.335(4) 


.176(4) 


.333(8) 



TABLE II: SU(2) + PhN model on V = 4 4 : Crossover parameters 



p 


rep/dist 




A 


P 


\[Aplv 






[A m , Am] 




p=.02 


.03 


.04 


.05 


.06 


.08 


.10 







F/SED 




.00933(1) 


.079(1) 


.117(1) 


.154(1) 


.194(1) 


.235(2) 


.311(2) 


.390(1) 


.376(0) 




F/LSD 


[.10, .30] 


.00929(0) 


.078(0) 


.117(0) 


.156(1) 


.194(1) 


.233(1) 


.311(2) 


.387(3) 


.375(0) 




A/SED 




.00617(0) 


.094(2) 


.145(2) 


.193(3) 


.244(4) 


.288(5) 


.381(8) 


.481(7) 


.377(2) 




A/LSD 


[.10, .50] 


.00618(0) 


.095(3) 


.142(4) 


.191(5) 


.237(6) 


.285(8) 


.383(15) 


.471(18) 


.374(4) 


0.5 


F/SED 




.01020(1) 


.071(1) 


.106(1) 


.143(1) 


.178(1) 


.213(1) 


.286(1) 


.353(2) 


.360(0) 




F/LSD 


[.10, .26] 


.01009(0) 


.072(0) 


.108(0) 


.143(1) 


.179(1) 


.215(1) 


.286(1) 


.359(4) 


.361(0) 




A/SED 




.00626(2) 


.094(2) 


.140(2) 


.185(3) 


.230(4) 


.274(4) 


.367(5) 


.446(10) 


.365(2) 




A/LSD 


[.10, .50] 


.00622(0) 


.094(3) 


.142(4) 


.189(5) 


.235(6) 


.283(8) 


.377(10) 


.468(28) 


.372(4) 


1.0 


F/SED 




.01150(5) 


.062(1) 


.094(1) 


.125(1) 


.156(1) 


.188(1) 


.255(2) 


.316(2) 


.338(0) 




F/LSD 


[.10, .30] 


.01120(0) 


.065(0) 


.097(0) 


.130(0) 


.162(1) 


.194(1) 


.259(2) 


.322(3) 


.343(0) 




A/SED 




.00631(3) 


.092(2) 


.140(2) 


.188(3) 


.233(4) 


.286(5) 


.376(7) 


.490(12) 


.373(2) 




A/LSD 


[.10, .50] 


.00632(0) 


.094(3) 


.138(4) 


.186(5) 


.231(6) 


.280(8) 


.374(14) 


.459(27) 


.369(4) 



TABLE III: x 2 /d-o.f. for fitting spectral distributions 



dist 


rep 


SU(2)+TBC 


SU(2)+PhN 


0=0 


0.5 


1.0 


1.5 


13=0 


0.5 


1.0 


SED 


F 


0.79-1.22 


0.54-1.00 


0.64-1.80 


0.86-3.13 


0.58-1.33 


0.63-1.50 


0.73-1.44 


A 


0.61-1.11 


0.36~1.24 


0.80-1.59 


0.79-1.31 


0.41-1.42 


0.65-1.88 


0.60-1.49 


LSD 


p* 


0.65-1.07 


1.11—1.53 


0.63-1.61 


0.89-1.41 


0.61-1.54 


0.55-1.29 


0.78-1.71 


A 


0.97-1.55 


0.63-1.09 


0.57-1.15 


0.60-1.17 


0.61-1.40 


0.80-1.22 


0.72-1.21 



* excludes F/LSD for TBC at p < 0.1, for which xV d -°- f =3 ~ 5. 
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FIG. 6: Smallest eigenvalue distributions of SU(2)+TBC L> KS in fund, (left) and adj. (right), at /3 = and tp = 0.01 ~ .06. 




FIG. 7: Microscopic level densities of SU(2)+TBC D KS in fund, (left) and adj. (right), at /3 = and <p = 0.01 ~ .06. 




FIG. 8: Level spacing distributions of SU(2)+TBC D KS in fund, (left) and adj. (right), at /3 = and ip = 0.01 ~ .06. 
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FIG. 9: Smallest eigenvalue distributions of SU(2)+PhN D KS in fund, (left) and adj. (right), at /3 = and p = 0.02 ~ .10. 




FIG. 10: Microscopic level densities of SU(2)+PhN D KS in fund, (left) and adj. (right), at /3 = and p = 0.02 ~ .10. 




FIG. 11: Level spacing distributions of SU(2)+PhN D KS in fund, (left) and adj. (right), at /3 = and ip = 0.02 ~ .10. 
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IV. LOW-ENERGY CONSTANTS 



A. chiral Lagrangian 



Below we briefly review the chiral Lagrangian description of QCD-like theories and sketch the identification of the 
low-energy constants with the parameters of random matrices. The effective low-energy Lagrangian for QCD-like 
theories with Np flavors of quarks in (pseudo)real representation, at finite chemical potential p and bare quark mass 
m is unambiguously fixed by the global symmetry alone (provided that p is much smaller than the p meson mass) 
and takes the form containing two phenomenological free parameters F and £ [B], 

C cS (Q) = ^F 2 trd^d^Q + 2F 2 ptrBQ^d Q - F 2 p 2 tr{BQ' i BQ + BB) - ^SmRetrMQ. (17) 

Here Q(x) is an SU(2Ai?) matrix- valued Nambu-Goldstone field, B = 03 ® liv F) M = 01 ® 1n f {i&2 ® 1jv>) f° r 
quarks in a real (pseudoreal) representation. F is the 'pion' decay constant and £ = Otptp) /Nf the chiral condensate, 
both measured in the chiral and zero-chemical potential limit m, p — > 0. If the theory is in a finite volume V = L 4 
and Thouless energy defined as E c ~ F 2 /EL 2 is much larger than m, the path integral is dominated by the zero-mode 
integration 

Z= I dQ exp[Vp 2 F 2 tT(BQ^BQ + BB) + -VZmRetrMQ) , (18) 
and the theory is said to be in the £-regime. In order to extract Dirac spectrum, one introduces fictitious bosonic 



quarks as well as fermionic quarks in the fundamental theory, leading to the graded group version of ( 18 ) on the 
effective theory side. For the actual computation one needs to parametrize the graded matrix Q in terms of its 
eigenvalues. Comparing the resulting expression (after analytic continuation n —¥ ifi and m — > iX) with the random 
matrix results ([6| and ([7]), the coefficients of chemical-potential and 'mass' terms in the exponents on both sides are 
readily identified as 4VF 2 fi 2 — 2tt 2 p 2 and VEX = irx. The latter is merely the definition of unfolded eigenvalues 
x = X/A, when combined with Banks-Casher relation £ = ir/AV that determines one of the low-energy constants S 
in terms of the mean spacing A of small Dirac eigenvalues. Eliminating the volume in favor of the level spacing, the 
former equation reads 



2 F 2 

VAp=\l-—n, (19) 

where the left hand side is a volume-independent combination. Accordingly one can determine another low-energy 
constant F 2 /£ from the slope of p,-\[Ap plots, preferably on lattices of various sizes. Note that in the parameter 



region V£|m| 3> 1, eq.(18) should approach the a model of nonchiral parametric random matrix ensembles [25], but 
we see no reason that the 'pion decay' or diffusion constant multiplying tr BQ^ BQ is affected. Accordingly, if the 
mean level spacing is approximately constant in a window in the very vicinity of the origin, one can determine F 2 /£ 
from the bulk correlation (namely the LSD) in that window. 
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FIG. 13: Low-energy constants F 2 and E for SU(2) quenched lattice gauge theory. Values are multiplied by 2/3 for the 
adjoint. 

B. SU(2) + twisted BC 

As the imaginary chemical potential is equal to the flux per link, p = (2tt /L)(p, F 2 /E can be extracted from the 
ip-\/~Ap plots for the SU(2)+TBC model. In Table I we exhibit the values of F 2 /E determined from the slopes of 
these tp-yAp plots, with all the numerals in the lattice unit. We also list the mean level spacings (A(0) for the 
SED rows and Zl([A m , Am]) for the LSD rows) and the chiral condensate E. The pion decay constant F 2 is then 
obtained by multiplying F 2 /E and E in each row. Coupling dependence of these low-energy constants as determined 
from respective spectral distributions are summarized in Fig. 13. In order to offset the number of components in the 
gauge multiplet, the constants for adjoint quarks are multiplied by 2/3 in the plot. We notice that except for the 
SU(2)-fundamental at j3 = 0.5 and 1.5, two values of F 2 , one determined solely from the zero-virtuality correlations 
(SED) and the other from the bulk correlations (LSD), agree within the one-cr error bars. Although it is admittedly 
difficult to calibrate the systematic deviations involved in employing the bulk correlation for the determination of F 
(part of which could originate from non-zero modes of the effective theory or from the nonuniformity of the mean level 
spacing within the eigenvalue window (see Fig. 5, left) due to the smallness of the lattice), these numerical agreements 
could justify the latter method using LSDs a posteriori. 

At /3 = 0, both low-energy constants agree between fundamental and adjoint representations. This observation is 
consistent with the fact that in the strong coupling limit the low-energy constants are identical to the Sp(-/V) and 
O(A^) lattice gauge theories at large N [51] . which share the same antiunitary symmetries as SU(2)-fundamental and 
adjoint. We also notice that due to the onset of chiral symmetry restoration at f3 = 1.5, the microscopic level density 
Ri(x) for SU(2)-fundamental starts to deviate from the random matrix result at x > 2 (Fig. 14, left). Thus the range 
of x available for fitting the SED is rather limited, especially on a lattice as small as V — 4 4 . On the other hand, the 
LSDs near the origin, not being directly sensitive to the chiral symmetry, arc still httable to the random matrix result 
without noticeable deviation throughout the plotted region (Fig. 14, right), and the values of V Ap determined from 
the LSD retains the linear dependence on (p. Therefore we consider the fitting of LSDs to the parametric nonchiral 
random matrices to be an efficient method of extracting F 2 from small lattices. 

C. SU(2) + phase noise 

The strength of phase randomness p in the SU(2)+PhN model is not related to the imaginary chemical potential, 
and thus cannot be directly used to determine F 2 from the response of the Dirac spectra. Thus we merely list in 
Table II the slopes of p-\f~Ap plots as a counterpart of p-\[A~p plots. We however noticed, rather unexpectedly, that 
for the SU(2) theory at the strong coupling limit f3 = the 'conversion ratio' between the phase randomness p and 
the imaginary chemical potential p is unity within numerical error (see the v Apj ' p column in Table I and the VAp/p 
column in Table II). This fact will have to be accounted for analytically. They gradually disagree for increasing j3. 

Finally, we need to check the volume dependence of the crossover parameter, which should scale as p = 
(V2F p/ir)\/V. Only for the purpose of varying the volume of the lattice in small steps, we adopted the two- 
dimensional toy model of SU(2)+PhN at P = 0. N con f — 22000 independent configurations are generated for each set 
of parameters. First we confirmed that the model shares the features in four-dimensions, i.e. all spectral distributions 
fit nicely to the parametric random matrices, and V~Ap scales with the randomness p on a lattice of fixed size V = 16 2 
(Fig. 15). The eigenvalue windows for sampling levels spacings are chosen to be [A m ,AM] = [0.30,0.65] (fundamental) 
and [0.16, 0.60] (adjoint). Then the p parameter is measured from the level spacing distributions on lattices of size 



■ — 1 


. Fnd/SED 




. Fnd/LSD 




. Adj/SED 




. Adj/LSD 






. Fnd/SED 
. Adj/SED 



15 




FIG. 14: Microscopic level density (left) and level spacing distribution (right) of SU(2)+TBC D KS in fundamental represen- 
tation, at /3 — 1.5 and tp = 0.01 ~ .06. 



0.20 




FIG. 15: (left) Crossover parameter \fKp for SU(2)+PhN model at /3 = and p — 0.02 ~ .10. Lattice size: 16 2 , number of 
configurations: 22000. Eigenvalue windows: [A m , Am] = [0.30,0.65] (F), [0.16, 0.60] (A). Real lines are linear fits to the data. 

FIG. 16: (right) Volume dependence of the crossover parameter p for SU(2)+PhN model at f3 = and p = 0.03, determined 
from LSDs. Lattice size: 10 2 ~ 24 2 . Other parameters are common to Fig.15. 

V = 10 2 ~ 24 2 at fixed phase randomness p — 0.03. p is indeed seen to scale as expected, linearly increasing in \JV 
(Fig.16). 



V. CONCLUSIONS 

In this paper we have computed the level spacing and smallest eigenvalue distributions of Hermitian random matrices 
in the G(0,S)E-GUE and chG(0,S)E-chGUE crossover. The results are shown to be perfectly fittable to the adjoint 
and fundamental staggered Dirac spectra of quenched SU(2) lattice gauge theory whose antiunitary symmetry is 
weakly violated by twisted boundary condition or noisy phases. This leads to the precise determination of the pion 
decay constant F (in the chiral and zero-density limit) from the Dirac spectral data. This method, feasible on a 
small-size lattice, has an advantage over the conventional method of measuring the decay rate of axial correlators, 
which requires a large temporal dimension. 

Our treatment is complementary to the previous approach of determining F of two-color QCD from its Dirac 
spectrum: Akemann and collaborators have concentrated on the real chemical potential and measured the response of 
Dirac eigenvalues that permeate into the complex plane, and fitted the lattice data to the non-Hermitian parametric 
chiral random matrices 18] . On the other hand we considered the twisted boundary condition, i.e. the imaginary 
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chemical potential and measured the response of Dirac eigenvalues that crossover within the real axis, and fitted 
to the Hermitian parametric chiral random matrices. In terms of chiral Lagrangian the difference is solely in the 
sign of fj, 2 , and the integral formulae involving (original or modified) Bessel functions Q, Q originated from the a 
model action are shared by the two. Similar complementary treatments were applied for three-color QCD at real and 
imaginary isospin chemical potential, which correspond to non-Hermitian chiral random matrices |31j and Hermitian 
chiral random matrices in chGUE-chGUE crossover [32, 38, 52 respectively, resulting in a successful determination 
of F. Combined with the results reported in this paper filling the missing pieces, the fact that the Dirac spectral 
statistics' in all three cases (SU(2)-fund.+/x, SU(2)-adj.+/z, SU(3)-fund.+/ii so ) agree perfectly with the predictions 
from corresponding identical zero-mode-approximated chiral Lagrangians in both regions of signs of /j, 2 constitutes a 
solid evidence for the validity of analytic continuation in the /i-plane, which is much required for the actual physics 
i.e. three-color QCD at real baryon number chemical potential. 

We consider the use of imaginary chemical potential has a practical advantage for the following reason. In order 
to fit non-Hermitian Dirac spectra to the non-Hermitian random matrix result, one usually project the complex 
eigenvalues either to the real or imaginary axis, and this projection could blur the fitting. Freer two-dimensional 
motion of complex eigenvalues is likely to lead to large statistical fluctuation. On the other hand, confining the 
eigenvalues within the real axis averts such issues, yielding a precise fitting shown in Figs. 6^11. 

Sharpness of the fitting functions chosen, pi(s) and P(s), is also in our advantage. Note that the use of Wigner 
surmise (random 2x2 matrices) of pi(s) and P(s) for the parametric random matrices |53 [ 154 ) is an uncontrolled 
approximation and is not suited for the precise determination of F, especially from fitting in the range of s > 2 
[55], 5 The advantage of our treatment is that exponentially-fast, uniform convergence is guaranteed for Nystrom-type 
method at increasing order. 

Our next obvious step is to include weakly coupled QED in the simulation, rather than the phase noise treated in 
this paper. Our preliminary study shows that Dirac spectra of SU(2)xU(l) quenched gauge theory are again fitted 
well to the parametric random matrix predictions, and we are currently accumulating numerical data on lattices 
of larger size than the current paper. This two-color QCD + QED model could be of interest to the lattice gauge 
community in which the three-color QCD + QED simulation in pursuit of precise measurement of isospin-related 
observables has attracted attention recently 56 , although we are exploiting the very difference of symmetry of 
SU(2)xU(l) as compared to SU(3)xU(l). Other possible extensions are: 

• In order to reduce the statistical error of fitting SEDs p\(s) for which the spectral averaging is not applic able , 
one could simultaneously use the k th smallest eigenvalues distribution pk{s),k — 2,3,... (computable by eq.(ll)) 
alongside with pi(s) for fitting p, and look for the overlapping of their error bars. 

• Extension of our treatment to Dirac spectrum in a nontrivial SU(2) gauge field topology is straightforward: on the 
lattice side one should measure the overlap Dirac spectrum |52) . whereas on the random matrix side the indices of 
the Bessel J functions in eqs.([6| and ^ are to be incremented by the topological charge v. 

• Introduction of dynamical quarks is interesting beyond the obvious reason of approaching more 'realistic' models 
of QCD; the weakly symmetry-violating U(l) gauge randomness couples to the SU(2) gauge randomness through 
the quark loops (fermion determinant) only in that case. This correlation between perturbed and perturbing can 
possibly bring nontrivial distortion to the relationship between the bare U(l) coupling constant and the crossover 
parameter, i.e. the pion decay constant. Introduction of finite quark masses on the random matrix side could 
become cumbersome, but in the light of individual eigenvalue distributions for chGUE-chGUE crossover (the latter 
of Ref.[32j, the computation is still feasible. 

• Our result suggests a possibility of an exotic continuum limit, in which the rate of decoupling of U(l) gauge field is 
adjusted to the rate of approaching the thermodynamic limit, while keeping the p 2 F 2 fixed. Albeit a rather artificial 
limit, this procedure might define a theory in which a nontrivial effect is induced to the interaction of Nambu-Goldstone 
bosons by the decoupling gauge interaction through the violation of antiunitary symmetry. The asymptotic slavery of 
U(l) is not an essential obstacle to this possibility, because any gauge field in the complex representation, namely the 
fundamental of SU(./V' > 3), is equally suited as a pertinent perturbation violating the antiunitary symmetry of SU(2). 

In forthcoming papers we wish to complete the project started here by covering up the directions listed above. 



5 Take for instance the asymptotic behavior of the level spacing distribution of GOE, log P(s) ~ —7s 2 (s 2> 1). The Wigner-surmised 
value 7 = 7r/4 is larger than the exact value (ir/4) 2 by 0.17, which brings a fatal 100% error at s ~ 2. 
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